Nonergodic Brownian Dynamics and the Fluctuation-Dissipation Theorem 
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Nonergodic Brownian motion is elucidated within the framework of the generalized Langevin equa- 
tion. For thermal noise yielding either a vanishing or a divergent zero- frequency friction strength, the 
non-Markovian Browninan dynamics exhibits a riveting, anomalous diffusion behavior being charac- 
terized by a ballistic or possibly also a localized dynamics. As a consequence, such tailored thermal 
noise may cause a net acceleration of directed transport in a rocking Brownian motor. Two no- 
table conditions for the thermal noise are identified in order to guarantee the fluctuation-dissipation 
theorem of first kind. 
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The phenomenon of Brownian motion has assumed 
a fundamental and influential role in the development 
of thermodynamical and statistical theories and contin- 
ues to do so as an inspiring source for active research 
in various fields of natural sciences 0. The Brown- 
ian dynamics can conveniently be described by a gen- 
eralized Langevin equation (GLE), reading: mv(t) + 
to Jo* 7(t, t')v(t')dt' + d x U{x, t) = e(t). The thermal ran- 
dom force e(t) is assumed to be uncorrelated with the ini- 
tial velocity. The memory friction 7(t, t') is typically re- 
lated to the correlation of the random forces |2|,|3|. Kiibo 
has addressed the common behavior of a classical equi- 
librium bath by setting (i) (e(t)s(t')) — mfesT7(£, £/) 
with (ii) 7(£, t') = j(\t — t'\) being time-homogeneous. 
The one-sided Fourier transform consequently then obeys 
iZepy(u;)] > 0, for real-valued v. Here fee is the Boltz- 
mann constant and T denotes the bath temperature. The 
validity of the GLE is thus restricted to the case of a sta- 
tionary thermal noise obeying the above two conditions; 
e.g. see in Refs. 0, S 01 ■ The complex- valued mo- 
bility £b(u>) and the memory-friction j(u>), respectively, 
are given in terms of the Fourier-Laplace transforms of 
the time-homogeneous correlations of the particle veloc- 
ity v(t) and the environmental noise e(t), respectively, 
as 
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fi{uj) = —— dtexp(-iu]t)(v(t )v(t +t)), (1) 
kb± Jo 

1 f°° 

»»7(w) = r-= / dtexp(-iut)(e(t )e(t +t)). (2) 
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To differentiate between these two relations (1) and (2), 
one refers to the first relation as the FDT of the first kind 
and the second one as the FDT of the second kind, re- 
spectively |2j, |3|] • The former characterizes the relaxation 
of the average particle velocity, while the latter yields the 
generalized susceptibility, describing the response of the 
bath degrees of freedom. 

These two FDT relations play a central role in the 
theory of thermal Brownian motion; its validation and 



generalization presents a timely subject that is presently 
hotly debated, both within theory and experiment. 

With this work we aim at demonstrating that Kubo's 
conditions for the equilibrium bath are generally not 
complete within the framework of linear response the- 
ory. This is so, because the existence of anomalous dif- 
fusion has not been considered in the original treatment 
by Kubo and others. We thus not only seek to complete 
set of conditions for general Brownian motion to obey an 
ergodic behavior but, more importantly, we like to dwell 
on rich behavior that emerges if ergodicity does not hold. 
Throughout the following we use the nomenclature of er- 
godicity in the context of the theory of stochastic pro- 
cesses, i.e. for open systems that are subjected to noise: 
Then the process is ergodic if the temporal long-time av- 
erage of a state function equals its stationary ensemble 
average 5J. Thus, we work beyond the well-known case 
of isolated systems obeying an area preserving, Hamil- 
tonian dynamics Put differently, our main objec- 

tive is to extend the theory of classical Brownian motion 
by focussing on the intricacies of a pos sible non-ergodic, 
non-Markovian dynamics H H [H GH E2 In par- 
ticular, the asymptotic long-time statistical probability 
density will be termed ergodic if it approaches a station- 
ary state that is independent of the choice of the initial 
preparation. 

Noisy dynamics exhibiting nonergodic behavior. As 
we shall show below, systems possessing via the FDT 
of the second kind no finite friction at zero frequency 
typically exhibit an non-ergodic behavior. Thus, there 
is an abundance of open system dynamics for which er- 
godicity is not granted. The thermal fluctuations result 
typically from a coupling to environmental degrees of 
freedom. In the cases of bi-linear system-bath coupling 
these are characterized by a corresponding spectral den- 
sity of bath modes, J(u>), obeying i2epy(a;)] = J(u))/mu> 
E3. E UM- Thus, if the coupling to low frequency 
modes is weak, as it is the case with optical-like phonons 
13 ! or broadband colored noise 01 j or a l so f° r the 
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celebrated case with a blackbody radiation field of the 
Rayleigh-Jeans type 01 the static friction vanishes at 
zero frequency u> = 0. Then, no efficient mechanism 
necessary for the compliance of the ergodic behavior is 
at work. Other situations that come to mind involve 
the vortex diffusion in magnetic fields |l9j |. or diverse 
other open dynamics with an inherent velocity-dependent 
system-bath coupling |l8|, l20i l2l| . Note also that a typi- 
cal solid state Drude bath spectrum is proportional to uj 3 , 
thus yielding as well a vanishing, zero- frequency friction 

DIE}. 

Conditions for the validity of the FDT. The GLE was 
originally derived by Mori by use of the Gram-Schmidt 
procedure; it was re-obtained by Lee using the recur- 
rence relations method j^. Starting out from the well- 
known system-plus-oscillator-reservoir model detailed in 
Ref. |2j| , the equations of motions directly yield the GLE 
with a stationary, correlated thermal noise e(i). The ini- 
tial coordinate and velocity of each oscillator are assumed 
to be distributed according to thermal equilibrium. We 
now add two additional mixing conditions to the above 
given, well-known relations (i) and (ii). These read: (iii) 
lim s ^o[ s _1 7( s )] ~^ 00 an d (iv) lim s ^o[s7(s)] = 0, where 
7(s) is the Laplace transform of the memory friction ker- 
nel that enters the GLE. Our terminology will be as fol- 
lows: If one of the above given four conditions (i)-(iv) is 
not satisfied we term the bath - in the specified order (i) 
- (iv): (i) a no neq uilibrium bath p^ . (ii) a relaxing (or 
ageing) one 0,|25[, (hi) a ballistic one 0,0], or (i y ) a 
nonergodic-bath, respectively. The notion embracing the 
latter two cases (iii-iv) will also be referred to as weak 
nonequilibrium heat baths. 

For the force-free GLE, i.e. with U(x) — 0, the two- 
time velocity correlation function (VCF) reads, 



C vv (h,t 2 ) = {v 2 {0)}b 2 + 
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where r = \ti — ia| and the explicit forms of the two relax- 
ation functions A and S arc detailed in Ref. 20] . Note 
that depending on the choice of initial preparation this 
correlation generally is not time-homogeneous. Herein, 
we indicate by {■ ■ •} the average with respect to the ini- 
tial preparation of the state variables. The relevant, gen- 
erally non-vanishing quantity b is given by: 



&=[l + limft(a)/«)]- 

s—>0 



(4) 



The non-equilibrium, time-inhomogeneous relaxation 
part A(ii,£a) obeys A(t\ — » oo,£2 — » oo) = 0. The 
asymptotic equal-time dynamics is found to read |20| : 
{(v 2 )} st = {v 2 (0)}b 2 + keTm- 1 ^ - b 2 ). Moreover, 
the asymptotic velocity correlation function emerges as 
C vv st = bk B T/m for the case that we set the initial 
velocity variance in accordance with the thermal equi- 
librium value {v 2 (0)} = kBT/m. These results follow 



for b ^ 0, thus implying that the condition (iii) is not 
satisfied (case with a ballistic bath). This in turn re- 
quires that 7(0) = J °° 7(i)di = 0, meaning a vanishing 
effective friction at zero frequency, or equivalently a van- 
ishing of the spectral density of the noise e(t) at u> = 0, 
see in (2) Tnis causes a breakdown of the FDT of 

the first kind because of the preparation-dependent, i.e. 
i>(0)-dependent asymptotic result. This is so although 
the FDT of the second kind in (2) is valid. 

Next, we discuss the situation with a nonergodic bath 
with lim s ^o[s7(s)] =/= 0. Towards this goal, we consider 
the solution v(t) of the GLE for a case where b ^ 
constitutes a stationary noise source that drives the force- 
free Brownian GLE-dynamics. Put differently, we set for 
the acting thermal (tailored) noise source 



£i(i) = cv(t) — c 



{v(0)}R(t) 



1 



R(t - t')e(t')dt' 



(5) 

where the parameter c denotes a coupling strength. The 
response function is R{t) = b + J2j res [R( s j)] exp(sji), 
with its Laplace transform reading R(s) = (s + 7(s)) . 
The set {sj} denote the nonzero roots of the equa- 
tion s + 7(s) = 0, and the probability density for u(0) 
is chosen as a Gaussian with zero-mean and variance 
{k 2 (0)} = k B T/m. In this case, the resulting memory- 
friction kernel corresponding to the effective thermal 
noise e^t) is given by 7 i(t) = c 2 {(v(t)v(0))}/(k B T), 
yielding 71 (0) — > 00, but lim,,^o[ s 7i ( s )] = c 2 b = fi- 
nite. This finding implies that the Brownian particle with 
b 7^ is immersed in a nonergodic bath. 

Substituting (5) into the GLE, we obtain an additional 
quadric potential with the frequency c 2 b and a friction 
kernel given by 71 (t) — j(t) — c 2 b. This leads to 7(0) — » 
00 so that jS(0) = 0. Therefore, the diffusion constant 
vanishes, meaning that the motion of a force-free particle 
remains bounded. This phenomenon of localization at 
finite temperature may occur for a particle dynamics that 
is coupled to a set of frozen environmental oscillators 
corresponding to the sub-Ohmic bath density J(u>), being 
proportional to u s in the limit of 5 = [14| . A realization 
provides a harmonic chain on a Bethe lattice possessing 
no diffusion thus acting as a frozen bath in this context 

M 

The case with an inherent ballistic bath. We will 
concentrate on nonergodic Brownian motion where the 
Kubo FDT of the second kind is assumed to be sat- 
isfied, the condition (iii), however, is not met. A col- 
ored noise that induces ballistic diffusion is harmonic 
velocity noise (HVN) e(t) 27]. The HVN obeys the 
Langevin equations: y = e, i = —Te — Vl 2 y + 
where £(t) is Gaussian white noise of vanishing mean 
with (£(*)£(*')) = 2T]T 2 k B TS(t - t'), r) is the damp- 
ing coefficient, T and denote the damping and fre- 
quency parameters. The second moments of y(0) and 
e(0) obey {y 2 {0)} = rjTn- 2 k B T and {e 2 (0)} = nTk B T. 
The Laplace and Fourier transformations of the damp- 
ing kernel function read j(s) — i]Ts/(s 2 + Ts + fl 2 ) and 
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ReftM] = i]T 2 uj 2 /[(n 2 -lu 2 ) 2 + T 2 lu 2 }, respectively. The 
latter corresponds to the spectrum of HVN which now 
vanishes at zero-frequency. In presence of the FDT of 
the second kind with the conditions (i) and (ii) satisfied, 
the complex- valued mobility now emerges as 



AM 



m 



m(tt 2 + r)T) 
T(n 2 + nT) + ioj(n 2 + rjT - T 2 
(ft 2 + nT - uj 2 ) 2 + r 2 w 2 



(6) 



where b = (1 + jjTQ^ 2 )^ 1 and S(uj) — lim t _ KX) (l — 
exp(— iu>t))/(iu>). Therefore, A(0) — > oo, yielding a di- 
verging diffusion constant D — fi(Q)k B T — > oo. 

This free non-Markovian Brownian motion can also be 
recast as an embedded, higher-dimensional Markovian 
process. The Fokker-Planck equation (FPE) for the prob- 
ability density P{x,v,w,u,y,e;t) corresponding to the 
set of Markovian LEs via introducing the appropriate set 
of auxiliary- variables (w, U, y, e) obeys dtP — LppP = 0, 
where Lpp denotes the associated FPE operator: 



Lpp = —V-— 



d _! d 

V— TO W — 

OX ov 
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+ (Tw + nTv + fl 2 y + u) — 
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~n 2 (w -s)— +e— 
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+ (Ts + n 2 y)^+T 1 T 2 k B T °~ 



dw 2 



d 2 

-riTk B T— 2 . 



(7) 



The mean total energy of the thermal HVN-driven 
force-free particle reads {(E(t))} — \m{{v 2 {t))} = 
\m[{(v(t)) 2 } + {({v(t) - (v(t))) 2 )}}. Their stationary 
values are given by first two terms in the r.h.s of Eq. 
(3). In particular, the first part describes the rem- 
nant initial kinetic energy of the particle, which is dis- 
sipated partly by the present bath because of {(v(t — > 
oo))} = bv(0). This part vanishes for the common case 
with b = 0. The second part denotes the energy pro- 
vided from the heat bath and is independent of the 
initial particle velocity; it does not relax, however, to- 
wards equilibrium. The absorbed power of the particle 
is F abg = rjk B T[l + T?/(2r) + 2n 2 /(r 1 r)}- 1 , being smaller 
than the equilibrium value rjksT '. This implies that the 
Brownian particle is hindered in re-gaining the energy 
from a thermal noise that lacks a finite, zero-frequency 
spectral weight. Our numerical results are depicted in 
the figure 1, which are obtained from simulating a set 
of Markovian LEs equivalent to the FPE of (7). The 
numerics fully corroborate our theoretical findings. 

Non-steady transport in ratchet potential. A most 
intriguing situation refers to Brownian motors |28[ 
when driven by non-ergodic Brownian motion. Take 
the case of thermal HVN driving a Brownian par- 
ticle in a rocking ratchet potential: U(x,t) = 
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FIG. 1: (color online). Behavior of diverse mean energies (see 
text) for a nonergodic, force-free Brownian particle vs. time 
t. The parameters used are m = 1.0, ksT — 1.0, n — 0.2, 
T = 5.0, n = 1.0 x(0) = 0, and (i) v(0) = V2 and (ii) 1.0, 
respectively, for the total energies given by the solid lines and 
the remnant energies by the dashed lines, from top to bottom. 
The open circles is the adsorbed energy from the bath for all 
cases and also the total energy for the v(0) — 0.0 case. 



Uq [sin(27rx) + Ci sin(47ra:) + C2 sin(67rx)J + A(t)x with 

U Q = 0.461, ci = 0.245, and c 2 = 0.04 j^. Here A(t) is 
a square-wave periodic driving force that switches forth 
and back between A(t) = Aq when 2nt p < t < (2n+ l)t p 
and A(t) = -A when (2n + l)t p < t < 2{n + l)t p . We 
next study whether a stationary non-equilibrium current 
results that can be put to a constructive use in order to 
direct, separate or shuttle particles efficiently [2j|. We 
thus research the resulting, time- and ensemble-averaged 
current with the condition (iii) not obeyed. 

Figure 2 depicts the time-dependent, accelerated mean 
velocity [{(«(<))} = (2t p )- 1 f* +2tp dt'{(v(t'))}} for var- 
ious strengths of the driving force as obtained via the 
simulation of the Markovian LEs corresponding to the 
FPE of (7), being supplemented by the ratchet forcing 
—U'(x,t). A startling finding is that the averaged veloc- 
ity increases linearly with time. This directed accelera- 
tion is presented by the slope (see dashed lines in Fig. 2) 
of the average velocity. For a weak rocking (small ^4o) the 
phenomenon of directed motion involves the surmounting 
of barriers. In contrast, for strong rocking the averaged 
displacement is related to the mean square displacement 
through the modified Einstein relation, i.e., 



^2/^2(0) = k B T eS 



(8) 



The ballistic diffusion coefficient is K2 = lim^oo = 
{(x 2 '(t)}} a =o I '(2t 2 ) , the linear mobility ^(A) = 
0) follows from the nonlinear mobility /i2(A)) — 
lim t ^oo{(x(t))} / (Aot 2 ), and the effective temperature is 
T cS :—T + b 2 (m{v 2 {0)}/k B - T). The dissipative accel- 
eration of a particle of mass m subjected to a constant 
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FIG. 2: (color online). Accelerated, (time and ensemble)- 
averaged Brownian motor velocity in a rocking ratchet that is 
driven by HVN. The used parameters are m — 1.0, ksT = 0.5, 
■q = 5.0, T = 22.0, Q, 2 = 40.0, t p = 25.0, A = 2.0, 4.0, 6.0, 
10.0, and 15.0, from top to bottom. The particles undergo 
a finite acceleration a (times 100) vs. temperature T, being 
depicted with the inset for Aq — 10.0. 

force F yields a — Fb/m, with < 6 < 1. This finding 
is therefore intermediate between a purely Newton dy- 
namics (with 6=1) and an ordinary Langevin dynamics 
(with 6 = 0). Note also that the acceleration assumes a 
non-monotonic function vs. temperature. 



Resume. We have studied the nonergodic Brown- 
ian motion occurring in so termed weak nonequilibrium 
baths, where the fluctuation-dissipation theorem of the 
second kind still holds. The nonequilibrium results are 
directly related to the ergodicity breaking. The latter 
is manifested by the break-down of the FDT of the first 
kind. The nonergodic behavior is either due to a van- 
ishing or a divergent zero-frequency spectral density of 
the thermal noise. Our theory produces two limits of 
abnormal diffusions, namely ballistic diffusion and a lo- 
calized behavior. For ballistic diffusion the effective fric- 
tion vanishes at zero frequency while it assumes infinity 
for localization. In order to assure the ergodic behavior 
of an equilibrium bath the usual conditions (i) and (ii) 
must be completed by two additional conditions: (iii) 
lim s ^o[s _1 7(s)] 00 an d (iv) lim s _»o[s7(s)] = 0, where 
7(s) is the Laplace transform of the memory friction ker- 
nel. Yet another riveting result is that the correspond- 
ing Brownian dynamics for a rocking Brownian motor 
exhibits a distinct, accelerated, nonstationary velocity, 
rather than the constant drift which typifies the situa- 
tion with normal, Ohmic dissipation. 

We are also confident that our present results will 
crucially impact other quantities of thermodynamic and 
quantum origin. Thus, this field is open for future studies 
that in turn may reveal further surprising findings. 
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